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Resonant phonon depopulation terahertz quantum cascade lasers based on vertical and diagonal lasing tran- 
sitions are systematically compared using a well established ensemble Monte-Carlo approach. The analysis 
shows that for operating temperatures below 200 K, diagonal designs may offer superior temperature perfor- 
mance at lasing frequencies of about 3.5 THz and above; however, vertical structures are more advantageous 
for good temperature performance at lower frequencies. 
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Holding great promise as coherent sources in the 
terahertz and infrared regime, quantum cascade lasers 
(QCLs) have been subject to continuous optimization 
with respect to temperature performance and output 
power. Particularly, room temperature operation in the 
terahertz regime has been a long-standing goal. Progress 
has been made by introducing resonant phonon depop- 
ulation (RPD) designs?^ reducing the number of wells 
per device period^ and by applying copper metal-metal 
waveguides.— Optimizations with respect to the diago- 
nality of the lasing transition resulted in a 3.9 THz QCL 
operating at a record temperature of 186 Kj^ However, 
efforts in improving the temperature performance of a 
similar vertical 3.2 THz structure by introducing some 
degree of diagonality failed^ The goal of this paper is 
to clarify these apparently contradicting results and to 
investigate the role of diagonality for the optimization of 
RPD terahertz QCLs with respect to temperature oper- 
ation. 

We have designed and analyzed three-quantum-well 
terahertz QCLs with different degrees of diagonality 
and frequencies between 1.8 and 4.8 THz. For the self- 
consistent modeling of QCLs, advanced semiclassical en- 
semble Monte Carlo (EMC)£~— or full quantum trans- 
port methods^ are commonly used. Our EMC simula- 
tion tool has been specifically developed for the anal- 
ysis of QCLs ; 10 ' 11 self-consistently including all rele- 
vant scattering mechanisms, namely, electron-electron 
(e-e), electron- longitudinal optical (LO) and acoustic 
phonon, and electron-impurity scattering. Since inter- 
face roughness (IR) highly depends on the growth pro- 
cess and its exact parameters are difficult to measure, 
this scattering mechanism is included phenomenologi- 
cally using typical parameter values For the design 
of GaAs/Al .i5Gao.85As THz QCLs, barrier heights be- 
tween 135-150 meV are commonly usedi - — For our sim- 
ulations we assumed a barrier height of 165 meV, corre- 
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sponding to 72% conduction band offset^ This some- 
what larger value was utilized to avoid sharp anticross- 
ings of the lower laser level with weakly-bound upper 
states in the second downstream injector, arising in some 
of the investigated designs and leading to well-known 
simulation artifacts. 7 Effective electron masses of 0.067 
in the wells and 0.076 in the barriers are assumed. We 
verified that the observed trends are robust with respect 
to the assumed barrier heights and IR values. Our EMC 
simulation provides self-consistent results for the spectral 
gain in the terahertz structures, where the gain lincwidth 
is extracted from the above mentioned scattering rates 
based on lifetime broadening;— we note that nondiag- 
onal correlation effects can only be considered in fully 
quantum mechanical approaches.^ 

For our analysis, we have designed QCLs with different 
degrees of diagonality at various operating frequencies. 
The degree of transition diagonality between upper and 
lower laser states is quantified in terms of the oscillator 
strength f osc , normalized to the value for the correspond- 
ing vertical design /™ c rt . Thus, for the vertical structures 
we have /osc/Zosc* = 1) corresponding to a 0% diagonal 
design. A design with / OS c//o S c c rt = 0.7 is referred to as a 
30% diagonal structure, etc. In order to make the sim- 
ulation results comparable, a special effort was made to 
keep laser designs very similar. In particular, the upper 
and lower laser level anticrossing energies with injector 
states (for injection and extraction, respectively) were 
kept in the range 1.5-1.6 meV and 3.3-3.6 meV, respec- 
tively. If we assume the Alo.15Gao.85As barrier height 
of 135 meV, these injection/extraction anticrossing ener- 
gies become 1.95-2.0 meV for the injection anticrossing 
and 4.25-4.6 meV for the extraction anticrossing, which 
is in line with the injection/extraction anticrossing val- 
ues used in the current state-of-the-art devices^— The 
energy separations between the upper and lower injector 
states were kept in the range 37-39 meV in all laser de- 
signs. These parameters combined with the desired laser 
operating frequency and diagonality of the laser tran- 
sition uniquely define the layer sequence for the struc- 
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TABLE I. Overview of the designed THz QCLs. All layer thicknesses are given in angstrom, and bold numbers indicate barriers. 
The underlined wells are doped with a sheet density of 2.7xl0 10 cm -2 in their 55 A-wide middle region. 



Freq. 



0% diagonal 



30% diagonal 



50% diagonal 



1.8 THz 46/98/31/76/43/161 

2.3 THz 48/95/27/73/42/157 

2.8 THz 48/94/24/72/42/157 

3.2 THz 48/96/20/74/42/161 

4.1 THz 47/99/15/73/40/164 

4.8THz 49/98/12/71/39/163 



70% diagonal 



49/94/35/78/47/161 
51/90/31/77/46/158 
51/89/28/78/46/159 
51/90/24/81/46/163 
51/88/19/83/44/164 
52/86/15/84/42/164 



51/92/39/80/48/161 
52/89/35/80/47/159 
52/87/32/80/47/159 
52/88/29/84/47/163 
52/86/23/87/45/164 
53/83/20/89/44/164 



52/91/48/80/49/161 
52/88/42/81/48/160 
53/86/39/81/48/159 
52/87/36/86/48/163 
53/85/32/89/47/164 
54/82/29/92/45/164 
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FIG. 1. (Color online) Conduction band diagram for the 
3.2 THz QCL with (a) 0% diagonality and (b) 70% diago- 
nality. Upper and lower laser levels are marked by bold solid 
lines. 



ture^ The layer thicknesses of the designed structures 
are summarized in TableQ] Experimentally, the 0% diag- 
onal structures were tested for the 2.3 THz, 3.2 THz, and 
4.1 THz designs£±2 Also copper double-metal 3.2 THz 
QCLs with an active region based on a 30% diagonal 
transition were tested experimentally. The devices op- 
erated up to 174 K,- which is lower than the maximum 
operating temperature of 178 K for similar devices based 
on a vertical 3.2 THz design. 5 

As we go from 0% diagonal to 70% diagonal structures, 
the wavefunctions get more localized in a single well as 
shown in Fig. Q]for the 3.2 THz structure. From the basic 
perspectives of THz QCL design, the diagonal laser tran- 
sition is expected to help improve the electron injection 
efficiency into the upper laser state, suppress electron 
leakage from the upper laser state to the downstream in- 
jector, and reduce the nonradiative electron scattering 
rate from the upper laser state; however, the diagonal 
laser transition results in smaller transition dipole mo- 
ment compared to a vertical transition.— & One may ex- 
pect that THz QCL structures based on vertical design 
may provide higher gain at lower operating temperatures, 
when electron injection into the upper laser state is effi- 
cient and the LO phonon scattering of thermally excited 
electrons in the upper laser stated is suppressed. How- 
ever, diagonal transitions may have advantages at higher 
operating temperatures. Because the parameters of the 
electron transport (injection efficiencies and lifetimes in 
various laser states) arc difficult to estimate analytically, 
we use EMC simulations to compare the peak gain of de- 
vices with various diagonalities at different temperatures. 

Figure [5] shows the temperature dependent peak gain 
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FIG. 2. (Color online) Dependence of the peak gain on the 
active region lattice temperature for different operating fre- 
quencies and diagonalities. The shown results are for opti- 
mum bias, which only weakly depends on temperature and 
has for each structure been determined at 150 K as to yield 
maximum gain. 



for the designed structures in Table HI as obtained by 
our EMC simulations. As expected, for increased diago- 
nality, we observe a decreased gain at low temperatures, 
but also a smaller gain degradation with temperature. At 
laser lattice temperatures around 150-200 K, representing 
the currently relevant range for temperature performance 
optimizations, diagonality does not offer an advantage 
for the low frequency structures that we studied. Specifi- 
cally, in case of the 3.2 THz design we see no improvement 
of the temperature performance for the 30% diagonal as 
compared to the 0% diagonal structure, which is consis- 
tent with experimental observations £ For higher lasing 
frequencies, diagonal designs offer clear advantages. The 
reason is that when the energy spacing between upper 
and lower laser level increases towards the LO phonon en- 
ergy (36meV in GaAs), scattering of thermally excited 
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FIG. 3. (Color online) An and A/ for different degrees of 
diagonality for the (a) 2.3 THz QCL and (b) the 4.1 THz QCL 
at lattice temperatures of 150 and 200 K, respectively. The 
bias is chosen as to maximize the peak gain. 

electrons in the upper laser level becomes very strong 
even at modest temperatures^ resulting in reduced pop- 
ulation inversion and thus a strongly decreased peak gain. 
This detrimental mechanism is suppressed in diagonal de- 
signs, which for high operating temperatures maintain in- 
creased population difference outweighing their reduced 
oscillator strength. For low frequency structures, where 
LO phonon scattering of thermally excited electrons be- 
tween the laser levels is less pronounced, diagonal de- 
signs offer an advantage only at higher operating tem- 
peratures, where the thermal activation of LO phonons 
is stronger. Thus, for the currently reached operating 
temperatures*^ we find improvements in peak gain only 
for designs operating above approximately 3.5 THz. 

For additional details on our simulation results, we 
show in Fig. [3] the population inversion An between the 
upper and lower laser level and gain bandwidth A/ (full 
width at half maximum) at maximum gain, for 2.3 THz 
and 4.1 THz structures of different diagonalities. We note 
that the peak gain scales with / osc An/A/. Overall, the 
2.3 THz structures (Fig. |3ja)) exhibit a higher A/ in our 
simulations than the 4.1 THz designs (Fig. EJb)), which 
is mainly due to increased Coulomb scattering as well as 
IR scattering for these structures. Higher An and, to a 
lesser extent, smaller A/ of more diagonal designs over- 
compensate the reduction of f osc in case of the 4.1 THz 
structures, however not for the 2.3 THz QCLs. For both 
structures, a decrease of A/ with increasing diagonality 
is observed. The gain broadening is related to scatter- 
ing events involving the laser levels. For diagonal lasing 
transitions, all the contributions get reduced except for 
IR scattering, which may go up or down. Furthermore, 
we observe an increase of A/ with temperature, which is 
due to enhanced LO phonon scattering. 

The effect of diagonality d, is further investi- 
gated by introducing the relative quantities An re i = 
An(d)/An(0%) and A/ re i = A/(d)/A/(0%). In Fig. 1 
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FIG. 4. (Color online) Temperature dependence for (a) 30% 
diagonality and (b) 50% diagonality of the relative inversion 
An rc i and relative inverse linewidth A/~j . The peak gain 
improvement is oc An^iA/^ . 



we plot An rc i and A/ ro | 1 as a function of temperature for 
d = 30% ( Fig. IHa)) and d = 50% (Fig.^b)). Compar- 
ing Fig.|4ja) to (b), we find that An re i increases with tem- 
perature. This trend is most pronounced for 4.1 THz de- 
vices with 50% diagonality. Very little change of A/"* in 
the temperature range shown in Fig. [5] indicates that di- 
agonal structures at elevated temperatures mainly profit 
from improved inversion. Finally, Fig. 2] clearly shows 
that diagonal designs offer a much stronger improvement 
in inversion for the 4.1 THz QCLs than for the 2.3 THz 
structures. 

To investigate the robustness of our results, simula- 
tions assuming different IR scattering parameters as well 
as simulations with a fixed lasing transition linewidth (of 
1 THz) were performed (not shown). They all confirm 
the previously obtained trend of the spectral peak gains 
for the different structures. We find that for the case 
of a fixed linewidth, the gain improvement for diagonal 
designs is only somewhat reduced. 

In summary, we have designed three-quantum-well 
THz QCLs at various wavelengths featuring different de- 
grees of diagonality, and have analyzed the benefits of di- 
agonal laser transitions for high temperature operation. 
Between 150 and 250 K, we find that the main advantage 
of diagonal structures is the considerably increased inver- 
sion; additionally, the reduced gain bandwidth of diago- 
nal transitions is beneficial. These effects can outweigh 
the reduced oscillator strength and provide advantage 
in QCL temperature performance. For designs operat- 
ing above approximately 3.5 THz, we find that diagonal 
structures offer advantages at operating temperatures be- 
low 200 K. For lower frequencies the advantages offered 
by diagonal designs become relevant only at operating 
temperatures of 200-250 K or higher. These simulation 
results provide a basis for a further optimization of the 
temperature performance of terahertz QCLs. 
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